A comparison of software for analysis of rare and common short tandem repeat (STR) variation using human genome sequences from clinical and population-based samples.

Short tandem repeat (STR) variation is an often overlooked source of variation between genomes. STRs comprise about 3% of the human genome and are highly polymorphic. Some cause Mendelian disease, and others affect gene expression. Their contribution to common disease is not well-understood, but recent software tools designed to genotype STRs using short read sequencing data will help address this. Here, we compare software that genotypes common STRs and rarer STR expansions genome-wide, with the aim of applying them to population-scale genomes. By using the Genome-In-A-Bottle (GIAB) consortium and 1000 Genomes Project short-read sequencing data, we compare performance in terms of sequence length, depth, computing resources needed, genotyping accuracy and number of STRs genotyped. To ensure broad applicability of our findings, we also measure genotyping performance against a set of genomes from clinical samples with known STR expansions, and a set of STRs commonly used for forensic identification. We find that HipSTR, ExpansionHunter and GangSTR perform well in genotyping common STRs, including the CODIS 13 core STRs used for forensic analysis. GangSTR and ExpansionHunter outperform HipSTR for genotyping call rate and memory usage. ExpansionHunter denovo (EHdn), STRling and GangSTR outperformed STRetch for detecting expanded STRs, and EHdn and STRling used considerably less processor time compared to GangSTR. Analysis on shared genomic sequence data provided by the GIAB consortium allows future performance comparisons of new software approaches on a common set of data, facilitating comparisons and allowing researchers to choose the best software that fulfils their needs.

Detection of SARS-COV-2 variants and their proportions in wastewater samples using next-generation sequencing in Finland.

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants may have different characteristics, e.g., in transmission, mortality, and the effectiveness of vaccines, indicating the importance of variant detection at the population level. Wastewater-based surveillance of SARS-CoV-2 RNA fragments has been shown to be an effective way to monitor the COVID-19 pandemic at the population level. Wastewater is a complex sample matrix affected by environmental factors and PCR inhibitors, causing insufficient coverage in sequencing, for example. Subsequently, results where part of the genome does not have sufficient coverage are not uncommon. To identify variants and their proportions in wastewater over time, we utilized next-generation sequencing with the ARTIC Network's primer set and bioinformatics pipeline to evaluate the presence of variants in partial genome data. Based on the wastewater data from November 2021 to February 2022, the Delta variant was dominant until mid-December in Helsinki, Finland's capital, and thereafter in late December 2022 Omicron became the most common variant. At the same time, the Omicron variant of SARS-CoV-2 outcompeted the previous Delta variant in Finland in new COVID-19 cases. The SARS-CoV-2 variant findings from wastewater are in agreement with the variant information obtained from the patient samples when visually comparing trends in the sewerage network area. This indicates that the sequencing of wastewater is an effective way to monitor temporal and spatial trends of SARS-CoV-2 variants at the population level.

Quality Control of an Isogenic H/N/KRAS-Less Mouse Embryonic Fibroblast Cell Line Panel.

Isogenic H/N/KRAS-less mouse embryonic fibroblast (MEF) cell lines have been developed to assist in cell-based assays of RAS inhibitors. The quality control assessment of a panel of these isogenic MEFs is described here, with a focus on ensuring the proper insertion of the desired mutant RAS transgene, a determination of gene copy number, and an investigation of potential off-target mutations which could lead to phenotypes which are undesired in downstream experiments. Using this suite of quality control tools, a MEF cell line can be readily validated, and researchers can be assured of the rationale for an observed phenotype.

Case report: Primary vulvar adenocarcinoma of mammary gland type-its genetic characteristics by focused next-generation sequencing.

Mammary-like vulvar adenocarcinoma (MLVA) is an exceedingly rare subtype of vulvar adenocarcinoma that shares features with mammary gland tissue. Due to its rarity and lack of consensus, MLVA presents diagnostic challenges to pathologists. We present the case of a 59-year-old female with an ulcerated mass on the right side of the external genitalia, diagnosed as MLVA. Comprehensive immunohistochemistry (IHC) and gene sequencing studies were performed to characterize the tumor. IHC analysis revealed triple expression of hormonal receptors (estrogen receptor, progesterone receptor, and HER2), supporting the mammary gland origin of the tumor. Gene sequencing identified unique genetic mutations associated with the expression of hormonal markers. One fusion gene (ERBB2-NAGLU) has not been reported in any tumors, and other mutations with unique mutation types have not been previously reported in MLVA. Our findings shed light on the molecular characteristics of MLV and may help improve the diagnosis and treatment of this rare type of vulvar adenocarcinoma.

Genotyping Plasmodium falciparum gametocytes using amplicon deep sequencing.

BACKGROUND: Understanding the dynamics of gametocyte production in polyclonal Plasmodium falciparum infections requires a genotyping method that detects distinct gametocyte clones and estimates their relative frequencies. Here, a marker was identified and evaluated to genotype P. falciparum mature gametocytes using amplicon deep sequencing. METHODS: A data set of polymorphic regions of the P. falciparum genome was mined to identify a gametocyte genotyping marker. To assess marker resolution, the number of unique haplotypes in the marker region was estimated from 95 Malawian P. falciparum whole genome sequences. Specificity of the marker for detection of mature gametocytes was evaluated using reverse transcription-polymerase chain reaction of RNA extracted from NF54 mature gametocytes and rings from a non-gametocyte-producing strain of P. falciparum. Amplicon deep sequencing was performed on experimental mixtures of mature gametocytes from two distinct parasite clones, as well as gametocyte-positive P. falciparum field isolates to evaluate the quantitative ability and determine the limit of detection of the genotyping approach. RESULTS: A 400 bp region of the pfs230 gene was identified as a gametocyte genotyping marker. A larger number of unique haplotypes was observed at the pfs230 marker (34) compared to the sera-2 (18) and ama-1 (14) markers in field isolates from Malawi. RNA and DNA genotyping accurately estimated gametocyte and total parasite clone frequencies when evaluating agreement between expected and observed haplotype frequencies in gametocyte mixtures, with concordance correlation coefficients of 0.97 [95% CI: 0.92-0.99] and 0.92 [95% CI: 0.83-0.97], respectively. The detection limit of the genotyping method for male gametocytes was 0.41 pfmget transcripts/µl [95% CI: 0.28-0.72] and for female gametocytes was 1.98 ccp4 transcripts/µl [95% CI: 1.35-3.68]. CONCLUSIONS: A region of the pfs230 gene was identified as a marker to genotype P. falciparum gametocytes. Amplicon deep sequencing of this marker can be used to estimate the number and relative frequency of parasite clones among mature gametocytes within P. falciparum infections. This gametocyte genotyping marker will be an important tool for studies aimed at understanding dynamics of gametocyte production in polyclonal P. falciparum infections.

Automating the Illumina DNA library preparation kit for whole genome sequencing applications on the flowbot ONE liquid handler robot.

Whole-genome sequencing (WGS) is currently making its transition from research tool into routine (clinical) diagnostic practice. The workflow for WGS includes the highly labor-intensive library preparations (LP), one of the most critical steps in the WGS procedure. Here, we describe the automation of the LP on the flowbot ONE robot to minimize the risk of human error and reduce hands-on time (HOT). For this, the robot was equipped, programmed, and optimized to perform the Illumina DNA Prep automatically. Results obtained from 16 LP that were performed both manually and automatically showed comparable library DNA yields (median of 1.5-fold difference), similar assembly quality values, and 100% concordance on the final core genome multilocus sequence typing results. In addition, reproducibility of results was confirmed by re-processing eight of the 16 LPs using the automated workflow. With the automated workflow, the HOT was reduced to 25 min compared to the 125 min needed when performing eight LPs using the manual workflow. The turn-around time was 170 and 200 min for the automated and manual workflow, respectively. In summary, the automated workflow on the flowbot ONE generates consistent results in terms of reliability and reproducibility, while significantly reducing HOT as compared to manual LP.

Towards accurate indel calling for oncopanel sequencing through an international pipeline competition at precisionFDA.

Accurately calling indels with next-generation sequencing (NGS) data is critical for clinical application. The precisionFDA team collaborated with the U.S. Food and Drug Administration's (FDA's) National Center for Toxicological Research (NCTR) and successfully completed the NCTR Indel Calling from Oncopanel Sequencing Data Challenge, to evaluate the performance of indel calling pipelines. Top performers were selected based on precision, recall, and F1-score. The performance of many other pipelines was close to the top performers, which produced a top cluster of performers. The performance was significantly higher in high confidence regions and coding regions, and significantly lower in low complexity regions. Oncopanel capture and other issues may have occurred that affected the recall rate. Indels with higher variant allele frequency (VAF) may generally be called with higher confidence. Many of the indel calling pipelines had good performance. Some of them performed generally well across all three oncopanels, while others were better for a specific oncopanel. The performance of indel calling can further be improved by restricting the calls within high confidence intervals (HCIs) and coding regions, and by excluding low complexity regions (LCR) regions. Certain VAF cut-offs could be applied according to the applications.

A comprehensive benchmark of graph-based genetic variant genotyping algorithms on plant genomes for creating an accurate ensemble pipeline.

BACKGROUND: Although sequencing technologies have boosted the measurement of the genomic diversity of plant crops, it remains challenging to accurately genotype millions of genetic variants, especially structural variations, with only short reads. In recent years, many graph-based variation genotyping methods have been developed to address this issue and tested for human genomes. However, their performance in plant genomes remains largely elusive. Furthermore, pipelines integrating the advantages of current genotyping methods might be required, considering the different complexity of plant genomes. RESULTS: Here we comprehensively evaluate eight such genotypers in different scenarios in terms of variant type and size, sequencing parameters, genomic context, and complexity, as well as graph size, using both simulated and real data sets from representative plant genomes. Our evaluation reveals that there are still great challenges to applying existing methods to plants, such as excessive repeats and variants or high resource consumption. Therefore, we propose a pipeline called Ensemble Variant Genotyper (EVG) that can achieve better genotyping performance in almost all experimental scenarios and comparably higher genotyping recall and precision even using 5× reads. Furthermore, we demonstrate that EVG is more robust with an increasing number of graphed genomes, especially for insertions and deletions. CONCLUSIONS: Our study will provide new insights into the development and application of graph-based genotyping algorithms. We conclude that EVG provides an accurate, unbiased, and cost-effective way for genotyping both small and large variations and will be potentially used in population-scale genotyping for large, repetitive, and heterozygous plant genomes.

Improving HLA typing imputation accuracy and eplet identification with local next-generation sequencing training data.

Assessing donor/recipient HLA compatibility at the eplet level requires second field DNA typings but these are not always available. These can be estimated from lower-resolution data either manually or with computational tools currently relying, at best, on data containing typing ambiguities. We gathered NGS typing data from 61,393 individuals in 17 French laboratories, for loci A, B, and C (100% of typings), DRB1 and DQB1 (95.5%), DQA1 (39.6%), DRB3/4/5, DPB1, and DPA1 (10.5%). We developed HaploSFHI, a modified iterative maximum likelihood algorithm, to impute second field HLA typings from low- or intermediate-resolution ones. Compared with the reference tools HaploStats, HLA-EMMA, and HLA-Upgrade, HaploSFHI provided more accurate predictions across all loci on two French test sets and four European-independent test sets. Only HaploSFHI could impute DQA1, and solely HaploSFHI and HaploStats provided DRB3/4/5 imputations. The improved performance of HaploSFHI was due to our local and nonambiguous data. We provided explanations for the most common imputation errors and pinpointed the variability of a low number of low-resolution haplotypes. We thus provided guidance to select individuals for whom sequencing would optimize incompatibility assessment and cost-effectiveness of HLA typing, considering not only well-imputed second field typing(s) but also well-imputed eplets.

Clinical application of metagenomic next-generation sequencing for the diagnosis of suspected infection in adults: A cross-sectional study.

Metagenomic next-generation sequencing (mNGS) has become an available method for pathogen detection. The clinical application of mNGS requires further evaluation. We conducted a cross-sectional study of 104 patients with suspected infection between May 2019 and May 2021. The risk factors associated with infection were analyzed using univariate logistic analysis. The diagnostic performance of pathogens was compared between mNGS and conventional microbiological tests. About 104 patients were assigned into 3 groups: infected group (n = 69), noninfected group (n = 20), and unknown group (n = 15). With the composite reference standard (combined results of all microbiological tests, radiological testing results, and a summary of the hospital stay of the patient) as the gold standard, the sensitivity, specificity, positive predictive value, negative predictive value of mNGS was 84.9%, 50.0%, 88.6%, and 42.1%, respectively. Compared with conventional microbiological tests, mNGS could detect more pathogens and had obvious advantages in Mycobacterium tuberculosis, Aspergillus, and virus detection. Moreover, mNGS had distinct benefits in detecting mixed infections. Bacteria-fungi-virus mixed infections were the most common in patients with severe pneumonia. mNGS had a higher sensitivity than conventional microbiological tests, especially for M. tuberculosis, Aspergillus, viruses, and mixed infections. We suggest that mNGS should be used more frequently in the early diagnosis of pathogens in critically ill patients in the future.

Measuring, visualizing, and diagnosing reference bias with biastools.

Many bioinformatics methods seek to reduce reference bias, but no methods exist to comprehensively measure it. Biastools analyzes and categorizes instances of reference bias. It works in various scenarios: when the donor's variants are known and reads are simulated; when donor variants are known and reads are real; and when variants are unknown and reads are real. Using biastools, we observe that more inclusive graph genomes result in fewer biased sites. We find that end-to-end alignment reduces bias at indels relative to local aligners. Finally, we use biastools to characterize how T2T references improve large-scale bias.

Identifying novel disease genes and revealing the pathomechanism of monogenic diseases.

Diseases are caused by genetic and/or environmental factors. It is important to understand the pathomechanism of monogenic diseases that are caused only by genetic factors, especially prenatal- or childhood-onset diseases for pediatricians. Identifying "novel" disease genes and elucidating how genomic changes lead to human phenotypes would develop new therapeutic approaches for rare diseases for which no fundamental cure has yet been established. Genomic analysis has evolved along with the development of analytical techniques, from Sanger sequencing (first-generation sequencing) to techniques such as comparative genomic hybridization, massive parallel short-read sequencing (using a next-generation sequencer or second-generation sequencer) and long-read sequencing (using a next-next generation sequencer or third-generation sequencer). I have been researching human genetics using conventional and new technologies, together with my mentors and numerous collaborators, and have identified genes responsible for more than 60 diseases. Here, an overview of genomic analyses of monogenic diseases that aims to identify novel disease genes, and several examples using different approaches depending on the disease characteristics are presented.

Targeted phasing of 2-200 kilobase DNA fragments with a short-read sequencer and a single-tube linked-read library method.

In the human genome, heterozygous sites refer to genomic positions with a different allele or nucleotide variant on the maternal and paternal chromosomes. Resolving these allelic differences by chromosomal copy, also known as phasing, is achievable on a short-read sequencer when using a library preparation method that captures long-range genomic information. TELL-Seq is a library preparation that captures long-range genomic information with the aid of molecular identifiers (barcodes). The same barcode is used to tag the reads derived from the same long DNA fragment within a range of up to 200 kilobases (kb), generating linked-reads. This strategy can be used to phase an entire genome. Here, we introduce a TELL-Seq protocol developed for targeted applications, enabling the phasing of enriched loci of varying sizes, purity levels, and heterozygosity. To validate this protocol, we phased 2-200 kb loci enriched with different methods: CRISPR/Cas9-mediated excision coupled with pulse-field electrophoresis for the longest fragments, CRISPR/Cas9-mediated protection from exonuclease digestion for mid-size fragments, and long PCR for the shortest fragments. All selected loci have known clinical relevance: BRCA1, BRCA2, MLH1, MSH2, MSH6, APC, PMS2, SCN5A-SCN10A, and PKI3CA. Collectively, the analyses show that TELL-Seq can accurately phase 2-200 kb targets using a short-read sequencer.

Direct Determination of the Structure of Single Biopolymer Molecules Using Nanopore Sequencing.

This review highlights operational principles, features, and modern aspects of the development of third-generation sequencing technology of biopolymers focusing on the nucleic acids analysis, namely the nanopore sequencing system. Basics of the method and technical solutions used for its realization are considered, from the first works showing the possibility of creation of these systems to the easy-to-handle procedure developed by Oxford Nanopore Technologies company. Moreover, this review focuses on applications, which were developed and realized using equipment developed by the Oxford Nanopore Technologies, including assembly of whole genomes, methagenomics, direct analysis of the presence of modified bases.

EpiSegMix: a flexible distribution hidden Markov model with duration modeling for chromatin state discovery.

MOTIVATION: Automated chromatin segmentation based on ChIP-seq (chromatin immunoprecipitation followed by sequencing) data reveals insights into the epigenetic regulation of chromatin accessibility. Existing segmentation methods are constrained by simplifying modeling assumptions, which may have a negative impact on the segmentation quality. RESULTS: We introduce EpiSegMix, a novel segmentation method based on a hidden Markov model with flexible read count distribution types and state duration modeling, allowing for a more flexible modeling of both histone signals and segment lengths. In a comparison with existing tools, ChromHMM, Segway, and EpiCSeg, we show that EpiSegMix is more predictive of cell biology, such as gene expression. Its flexible framework enables it to fit an accurate probabilistic model, which has the potential to increase the biological interpretability of chromatin states. AVAILABILITY AND IMPLEMENTATION: Source code: https://gitlab.com/rahmannlab/episegmix.

scifi-ATAC-seq: massive-scale single-cell chromatin accessibility sequencing using combinatorial fluidic indexing.

Single-cell ATAC-seq has emerged as a powerful approach for revealing candidate cis-regulatory elements genome-wide at cell-type resolution. However, current single-cell methods suffer from limited throughput and high costs. Here, we present a novel technique called scifi-ATAC-seq, single-cell combinatorial fluidic indexing ATAC-sequencing, which combines a barcoded Tn5 pre-indexing step with droplet-based single-cell ATAC-seq using the 10X Genomics platform. With scifi-ATAC-seq, up to 200,000 nuclei across multiple samples can be indexed in a single emulsion reaction, representing an approximately 20-fold increase in throughput compared to the standard 10X Genomics workflow.

Literature Review of Pathogen Agnostic Molecular Testing of Clinical Specimens From Difficult-to-Diagnose Patients: Implications for Public Health.

To better identify emerging or reemerging pathogens in patients with difficult-to-diagnose infections, it is important to improve access to advanced molecular testing methods. This is particularly relevant for cases where conventional microbiologic testing has been unable to detect the pathogen and the patient's specimens test negative. To assess the availability and utility of such testing for human clinical specimens, a literature review of published biomedical literature was conducted. From a corpus of more than 4,000 articles, a set of 34 reports was reviewed in detail for data on where the testing was being performed, types of clinical specimens tested, pathogen agnostic techniques and methods used, and results in terms of potential pathogens identified. This review assessed the frequency of advanced molecular testing, such as metagenomic next generation sequencing that has been applied to clinical specimens for supporting clinicians in caring for difficult-to-diagnose patients. Specimen types tested were from cerebrospinal fluid, respiratory secretions, and other body tissues and fluids. Publications included case reports and series, and there were several that involved clinical trials, surveillance studies, research programs, or outbreak situations. Testing identified both known human pathogens (sometimes in new sites) and previously unknown human pathogens. During this review, there were no apparent coordinated efforts identified to develop regional or national reports on emerging or reemerging pathogens. Therefore, development of a coordinated sentinel surveillance system that applies advanced molecular methods to clinical specimens which are negative by conventional microbiological diagnostic testing would provide a foundation for systematic characterization of emerging and underdiagnosed pathogens and contribute to national biodefense strategy goals.